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The non-local vacuum condensates of QCD describe the distributions of quarks and gluons in 
the non-perturbative QCD vacuum. Physically, this means that vacuum quarks and gluons have 
nonzero mean-squared momentum, called virtuality. In this paper we study the quark virtuality 
which is given by the ratio of the local quark-gluon mixed vacuum condensate to the quark local 
vacuum condensate. The two vacuum condensates are obtained by solving Dyson- Schwinger Equa- 
tions of a fully dressed quark propagator with an effective gluon propagator. Using our calculated 
condensates, we obtain the virtuality of quarks in the QCD vacuum state. Our numerical predictions 
are consistent with other theoretical model calculations such as QCD sum rules, Lattice QCD and 
instanton models. 
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I. INTRODUCTION 

Quantum mechanics dictates that even "empty" 
space is not empty, but rather filled with quantum 
fluctuations of all possible kinds. In many contexts, 
such as in atomic physics, these vacuum fluctuations 
are subtle effects which can only be observed by pre- 
cision experiments. In other situations, especially 
when interactions of sufficient strength are involved, 
the vacuum fluctuations can be of substantial mag- 
nitude and even "condense" into a non- vanishing 
vacuum expectation value of some quantum fields, 
called vacuum condensates. These vacuum conden- 
sates can act as a medium^^l, which influences the 
properties of particles propagating through it. 

An important example of such a vacuum conden- 
sate is the Higgs vacuum expectation value, which is 
introduced in the Standard Model of particle physics 
to generate the masses of quarks, leptons, and the 
gauge bosons {W^ , Z^) of the weak interaction. 
The vacuum expectation value of the Higgs fleld, 
{(j)) — 2A6GeV , is uniquely determined in the Stan- 
dard Model. The quark and lepton masses differ 



from one another only due to the different strength 
of the coupling of each fermion to the Higgs fleld. At 
the same time, the quark masses also receive addi- 
tional contributions from the quark and gluon con- 
densates in the QCD vacuum. In fact, the contribu- 
tion of the QCD vacuum condensates to the masses 
for the three light quarks (u,d,s) considerably ex- 
ceed the mass believed to be generated by the Higgs 
field[2l. 

The non-vanishing value of chiral quark vac- 
uum condensates signals the spontaneous break- 
ing of chiral symmetry in QCD, and quantitatively 
it is related to the pseudo-Goldstone bosons mass 
spectrum [^1. Due to non-perturbative effects of 
QCD, the vacuum of QCD has a nontrivial struc- 
ture. The vacuum condensates are very important 
in the elucidation of the QCD structure and in de- 
scription of hadron properties. If the vacuum acts 
as a medium and influences the properties of funda- 
mental particles and their interactions, its properties 
can conceivably change. This idea has important im- 
plications in many aspects of physics. 

The non-perturbative vacuum of QCD is densely 
populated by long-wave fluctuations of quark and 
gluon flelds. The order parameters of this compli- 
cated state are characterized by the vacuum matrix 
elements of various singlet combinations of quark 
and gluon flelds, such as 

(0|:g-g:|0), {0 \: q[tgsa^,G';:,^]q -.1 0), 
(0 h 97myW7mY9 :| 0)' 
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^,,0), (0|:r''^G;i,G^,G^^:|0), •••,(!) 

which are called vacuum condensates of QCD, where 

q{x) is the quark field, GJJ,^ represents the gluon field 
strength tensor with a being color index (a = 1,2, 
• • •, 8), and can be expressed as 

G^^ix) = d^Al{x) - d.Al{x) + gsr'"^Al{x)Al{xU2) 

A" in expression (1) is the SU{?>) Gell- Mann ma- 
trix, /"^"^ represent the 5C/c(3) structure constants, 
and Qs in Eq.(2) is the coupling constant related 
to the so-called QCD running coupling constant ag 

by as{Q) = A^ is the gluon field, a^,^ = 

■^{lulv — Ivln) in Euclidean space with 7^ being a 
Dirac Matrix. 

In QCD by condensates we mean the vacuum 
mean values (0 | Oj | 0) of the local operators Oj {x) , 
which arise due to non-perturbative effects. The lat- 
ter point is very important and nc^ecls clarification. 
When determining vacuum condensates one implies 
the averaging only over non-perturbative fluctua- 
tions. If for some operators Oj the non-zero vacuum 
mean value appears also in the perturbative theory, 
it should not be taken into account in determination 
of the condensate. In other words, when determining 
condensates the perturbative vacuum mean values 
should be subtracted in calculation of the vacuum 
averages. 

Separation of perturbative and non-perturbative 

contributions to the quark propagator, Sq{x), has 
some arbitrariness. For the nonperturbative propa- 
gator, deflned in Sect. II, one makes an expansion 
of vaciuim expectation values involving antiquark- 
quark fields that would vanish in a perturbative vac- 
uum, the local vacuum condensates. The nonzero 
local quark vacuum condensate (0 |: q{0)q{0) :| 0) is 
responsible for the spontaneous breakdown of chiral 
symmetry. The nonzero local gluon vaciium con- 
densate (0 |: G°j,G°j, :| 0) defines the mass scale of 
hadrons through the trace anomaly W . 

The non-local vacuum condensates (0 |: q{x)q{0) :| 
0) describe the distribution of quarks in the non- 
perturbative vacuum . Physically, this means that 
vacuum quarks have a nonzero mean-squared mo- 
mentum called virtuality. Indeed, the quark average 
virtuality is connected with the vacuum expectation 
values'^'^'*' which will be discussed in the Sect. III. 

Studying the quark virtuality is of paramount 
importance for present day particle and nuclear 
physics, since it is not only related to the property 
of QCD vacuum states but also to quark vacuum 
condensates. In this work, we study the quark vir- 
tuality in the QCD vacuum by calculating quark and 



gluon vacuum condensates. The quark vacuum con- 
densates are obtained by solving Dyson-Schwinger 
Equations (DSEs)[^l of fully dressed quark prop- 
agators with an effective gluon propagator under 
the constraints of an Operator Product Expansion 
(OPE)[i°l. In Sect. II, we briefly introduce the DSEs 
for a fully dressed quark propagator and formulae of 
various quarks vacuum condensates. In Sect. Ill 
the quark virtuality is deflned; and the correspond- 
ing formulism of the virtuality is also derived in this 
section. In Sect. IV, our numerical results on the 
quark virtuality are presented. Our concluding re- 
marks of this study are given in Sect. V. 



II. DSES FOR QUARK PROPAGATOR 

To study quark virtuality, we need to know two 
quark vacuiim condensates, and quark gluon mixed 
vacuum condensates. Therefore, we begin with a 
study of the quark propagators, which determine 
various quark condensates and quark gluon mixed 
vacuum condensates under the OPE constraints. 
The quark propagator is deflned by 



S,{x) 



(0 I T[q-{x)q\0)] I 0) 



(3) 



where q°'{x) {q''{x)) is a quark fleld with color a 
(6), and T is the time-ordering operator. The 
fully dressed quark propagator in Eq.(3) can be 
decomposed into a perturbative part and a non- 
perturbative part. In other words, one can write 
the quark propagatorl^^'^^i as 

S,{x) = S^^{x) + S^''{x), (4) 

where, expanding in the quark mass m/. 



and 



1 7 • x 
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■_{0\:q{x)q{0) :| 0) 



+M0 |: q{x)rq{0) :| 0)] + 



(6) 



in configuration space, with a sum over color. For 
short distances, the Taylor expansion of the scalar 
part of S^^{x), (0 |: q{x)q{0) :\ 0), reads 

(0|:g(x)g(0):|0) = (0|:g(0)g(0):|0) 



-(0|:g(0)[iff,aG(0)](?(0):|0) + 



(7) 
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In Eq.(7) the local operators of the expansion are the 
local quark vacuum condensates, the quark-gluon 
mixed condensate, and so forth. 

An important observation is that the inverse 
quark propagator in momentum space can also be 
written in Euclidean space as 



Sj\p) = ip-Af{p^) + Bf{p% 



(8) 



which is renormalized at space - like point /i^ ac- 
cording to v4y(/i^) = 1 and Bf{fi'^) = mf{n'^), 
with mf{iJ?) being the current quark mass at re- 
normalization point /i^. The subscript f va. Af and 
Bf stands for quark flavor u,d and s. 

Except for the current quark mass and pertur- 
bative corrections, the functions [Af{p'^) — 1] and 
Bf(jP) are non-perturbative quantities which we re- 
fer to as the vector and scalar propagator conden- 
sates, respectively. The DSEs (in the Feynman 
gauge) satisfied hy Af and Bf then can be written 
as the set of coupled equations!^' ^^1. 

[Af{s)-l]s=^glJ^ s'ds' sin' xD{s,s') 

Vss' A f (,s') cos X 



s'A'j{s') + B'j{s') 



dx. 



(9) 



Bf{s) = ^glj s'ds' j sm'xD{s,s'\ 
Bf{s') 



s'A}{s') + B}{s') 



dx. 



(10) 



where s = p'^ and g-^D{s,s') — g^D{s + s' — 
2-\/ss' cosx) is the dressed gluon propagator. Now, 
our task is to solve this set of coupled equations, 
Eqs. (9,10), and get the solutions Af{s) and Bf{s). 

One can solve the two coupling integral equations, 
Eqs. (9,10), using an effective gluon propagator such 
as 



abi 



47ra(s) 



(11) 



where a{s) stands for quark-quark interaction which 
can be, for example, well approximated by 



(12) 



X in Eq. (12) is the strength of the interaction, 
and A is its range parameter. The first term of 



Eq.(12) simulates the infrared enhancement and con- 
finement, and the second term matches to the lead- 
ing log renormalization group results. The parame- 
ter e can be varied in the range 1.0 — 2.5. We take 
e to be 2.0 in the present calculations. The strength 
parameter x ^nd the parameter A are determined 
by fitting the solutions of DSEs to the pion decay 
constant [^^1, and they are listed in table 1. 

Table 1 . Values of the strength parameter x a-nd 
range parameters A of the quark-quark interaction 
used in our present calculations. 



Set no. 


Range A 


Strength x 


Set 1 


0.40 GeV2 


1.84 GeV 


Set 2 


0.20 GeV2 


1.65 GeV 


Set 3 


0.02 GeV2 


1.50 GeV 



The QCD scale parameter A and the value of d 
with the flavor number Nf = 3 are given by 

A = 0.2GeV, d= 12/(33 -2Nf) = 12/27 (13) 

The non-local quark vacuum condensate (0 |: 
q{x)q{0) :\ 0) is then given by the scalar part of 
Fourier transformed inverse quark propagator'^^'^^l. 



(0|:g(x)(z(0):|0) 
(-4iV.) f 



3 

'4^ 



(27r)4p2^2(p2)+52(p2) 



sds 



sA}{s) + Bj{s) 



(14) 



where the color number Nc = 3. Usi ng t he expan- 
sion of the Ji Bessel function, 2Ji{\/sx^)/Vsx'^ = 
1 — sa;^/8 -I- and the definitions of the condensates 
given in Eq(7), one finds that the quark condensate 
is 

(O|:,-(0WO):|0) = -jL£,.-j^^£M^15) 

while the local quark - gluon mixed vacuum conden- 
sate, (0 |: q{0)[igsaG{0)]q{0) :| 0) is 



(0 |: m[i9s(TG{0M0) :| 0) 
-A rdss'—""^^'^ 



sAj{s) + Bj{t 



(16) 



Note that the upper limit of the integrals over s for 
the DSEs, Eqs(9,10), is infinity, while for the conden- 
sates, Eqs(15,16), there is a finite limit, So- That is 
because the effective gluon propagator, Eqs(ll,12), 
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provides a natural cutoff for the DSEs integrals, 
while the range for the condensates is given by the 
renormalization point, /it. For light quarks perturba- 
tive QCD begins to dominate nonperturbative QCD 
at about 3-4 GeV, so a renormalization point of 
= 10 GeV^ is expected. Therefore, we use So = 10 
GeV^. This is explained in detail in Sect. IV. 

A different derivation of quark - gluon mixed vac- 
uum condensate explicitly used a form for the nonlo- 
cal quark condensate (0 |: q{x)q{0) :| 0) = g{x^){0 |: 
q{0)q{0) :| 0) (see Ref.[12]), and used a model for 
g{x) to derive an expression similar to that obtained 
in the model used in Ref.[14], but with very differ- 
ent results for the mixed quark-gluon condensate. In 
Sect. IV, we use the result for g{x) from Ref.[12] to 
estimate errors in our estimate for , defined in the 
following section. 

Eqs.(15) and (16) will produce our numerical pre- 
dictions of local quark and quark gluon mixed vac- 
uum condensates, which will be then used to es- 
timate the nonzero mean squared momentum of 
quarks in the non-perturbative QCD vacuum. 



III. NONZERO MEAN SQUARED 
MOMENTUM OF QUARK IN QCD VACUUM 

The quantities /{v) were introduced to represent 
nonlocal condensates^^^'^^l. Their explicit form com- 
pletely determines the coordinate dependence of the 
condensates, and describes the virtuality distribu- 
tion of quarks in the non-perturbative vacuum. Its 
n-th moment is proportional to the vacuum expec- 
tation value of the local operator with the covariant 
derivative squared to the n-th power^^^l; 



1 (0|:g(J^)"g:|0) 
r(n + 2) {0\:qq:\0) 



,(17) 



where the covariant derivative = — iggA^^ 
with = A1^X''/2 and A" is a SU{3) Gell-mann 
matrix. It is natural to suggest that the vacuum 
expectation values in the right-hand side of Eq. (17) 
should exist for any n. The two lowest moments 
{n = 0,n — 1) give the normalization condition (n = 

0) and the average vacuum virtuality of quarks {n = 

1) , A^, respectively 



and 



f,{iy)du = 1, (18) 



/■°° 1 (0 |: qD\ :| 0) _ A^ 

/ vfq{v)dv = - = (19 

/o 2 {Q\:qq:\0) 2 



To illustrate the definition of quark virtuality in the 
non-perturbative vacuum, A^, let us now consider 
the Taylor expansion of the simplest gauge invariant 
condensate 



(0 |: q{Q)E{Q,x-A)q{x) :| 0) = (0 |: q-(0)(z(x) :| 0) 
E-^Mi---a^M„(0|:9-i^'^^---i?^"g:|0) 



= {Q\:qq:\Q) + ^{Q\:qD^q:\Q) 



(20) 



where E = Pexp[i Af^{z)dz^^] is the path-ordered 
Schwinger phase factor (the integration is performed 
along the straight line ) required for gauge invariance 
and A^(z) = A;,(z)AV2. 

The quantity A^, defined in Eg p^ 



{0\:qD^q:\0) 

(0|:g-g:|0) ■ 



(21) 



was introduced for an expansion of the nonlocal 
quark condensate^^I, which can be interpreted as the 
average virtuality of the vacuum quarks. Note that 
the operator (0 |: qD^q :| 0) can be presented in a 
different form: 



(0 
(0 
(0 



qD\ :| 0) = (0 |: qOf^D^q :| 0) 
qD^^g^^D-'q -.10) 

qp pq :] 0) - {0\:qD^^D'' a ^,q -.10) 



(22) 



-mliO |: qq :| 0) - -(0 |: qlD^' , D-']a^,q :| 0) 



1 



-m^(0 |: g<7 :| 0) + -(0 |: qiigsG'^" a^,]q :| 0), 



where p = J^Dfi, and we have used the identity 



(23) 



the equation of motion 

pq{x) = -im^qix), (24) 
and the definition of the field strength tensor 

\D^,,D,\=-igG^,. (25) 

Thus, the "average virtuality" of the vacuum quarks 
(0 |: qD'^q :| 0) is directly related to the "average 
vacuum gluon field strength" (0 |: q[igsG'^''a^^]q :| 
0). In many papers^^^l one can find the notation 

(0 |: q[igsG^''a,,,]q :| 0) = ml{0 |: qq :| 0). (26) 
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Using Eqs. (21,22,26), one can write the vacuum 
quark virtuality as 



rrin 



Namely, 

2 ^ 1 (0 |: q[igsG^"a^^]q :| 0) 
« 2 (0|:gg:|0) 



(27) 



(28) 



For light quarks, the mass term is so small that 
it can be neglected. Finally, we arrive at 



1 (0 |: g[ig,G^-a^.]g :| 0) 

2 (0|:gg:|0) 



(29) 



Eq. (29) is starting point of our calculations on 
quark virtuality. As we see from Eq. (29), in order 
to get Ag we have to calculate quark-gluon mixed 
vacuum condensate (0 |: q[igsG^''a^^]q :| 0) and two 
quark vacuum condensate (0 | : : | 0) by the use of 
Eqs. (15) and (16). 



Our results for two quark local vacuum conden- 
sates arc consistent with the predictions by Gall- 
Mann-Oakes-Renner relation (GMOR)!^^^^"! , (m„ + 
'md.){0 |: qq :| 0) = -^m^/^, where to„ and nid 
are current quark masses with value of m„ + = 
9.7McV[20l, and = 140McV, ,U = 93MeV are 
the mass and decay constant of a pion, respectively. 
Substituting these values into the GMOR relation 
produces {0 \: qq :] 0) = -COOSTGeV^, which is rea- 
sonably consistent with the u,d quark condensates 
shown in Table 2, within errors. 

Our theoretical results in Table 2 and 3 are also 
consistent with the empirical values used widely in 
QCD sum rulesPi], with {0 \: qq :] 0)"''^ ~ -0.013 
GeV^, as we obtained with Set 1. We also obtained 
a result consistent with the predictions of Lattice 
calculations'^^!. 

Using Eq.(29) and our numerical predictions of 
the two quark local vacuum condensates, and the 
local quark-gluon mixed vacuum condensates given 
respectively in tables 2 and 3, the quark virtuality, 
the nonzero mean squared momentum of quarks in 
non-perturbative QCD vacuum state, is given by 



IV. NUMERICAL RESULTS AND ERROR 
ESTIMATE 

Using the solutions of DSEs of Eqs. (9,10), Af 
and Bf, with three different sets of the quark-quark 
interaction parameters given in Table 1, leads to our 
following theoretical predictions for the local two 
quark vacuum condensates, local quark-gluon mixed 
vacuum condensates via Eqs. (15,16). The corre- 
sponding theoretical results are listed in Table 2-3. 

Table 2. The local two quark vacuum condensates 

of QCD, (0 |: qq :| 0)^, / stands for quark flavor 
and n denotes renormalization point, /Lt^=10 GeV^. 



Set No. 


{0\:qq: \ O)"''' 


{0\:qq:\Or 


Set 1 


-0.013(GcV)-'' 


-0.071 (GcV)-"^ 


Sen; 2 


-0.007cS(GcVj'' 


-0.0()8(G(>V)-' 


Set 3 


-0.0027(GeV)3 


-0.065(GeV)3 



Table 3. The local quark-gluon mixed vacuum 
condensates, (0 |: q[igsaG]q :| 0);^, /i^=10 GeV^. 



Set No. 


(0 |: q[igsaG]q :| O)"-"^ 


{0\:q[ig,aG]q:\Oy 


Set 1 


-0.015(GeV)5 


-0.186(GeV)5 


Set 2 


-O.OlO(GcV)^ 


-0.189(GeV)5 


Set 3 


-0.0078(GeV)5 


-0.193(GeV)5 



_ l(0|:g(0)[^g.ap.G^.^]g(0):|0)„,d 



2 (0|:g(0)g(0) 
= 0.57GeV^ 



0)u,d 



(30) 



for u, d quark. Set 1. 
For s quark, we obtained 



1 (0 |: ?(0)[i5«<7^.G°,^]?(0) :| 0), 



(0|:g(0)g(0):|0), 



= l.SlGeV^, 



(31) 



for Set 1. 

In the summary of this section, our predictions of 
X"^ ^ and A^ with three different sets of quark-quark 
interaction parameters given in table 1 are listed in 
Table 4. 



Table 4. The virtualities of three light quarks u, d 
and s in the QCD vacuum state. 



Set No. 




A2[GeV]2 


Set 1 


0.57 


1.31 


Set 2 


0.68 


1.38 


S(;t 3 


1.43 


1.48 



All our theoretical results are in an acceptable 
range'^^l of A^ between 0.4 ~ 2.50Gey^. For exam- 
ple, for u and d quarks the standard QCD sum rule 
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estimationP4] gives A^^ = 0.4 ± O.lGeV^, the QCD 

sum rule analysis of pion form factor^^l produces 
A«_^ = OJOGeV^, and Lattice QCD calculations I^e] 

predicts A^^ = O.SSGeV^. For s quark, Lattice 
QCD P6] gives A2 = 2.50GeV'^, and the instan- 
ton model predictionl^^l^ = lAOGeV"^. The 
^ud ~ ^■^'•^ ^ larger than that for set 

1 (0.57) and set 2 (0.68). The reason is that the 
range parameter A is an order of magnitude smaller. 
Therefore we observe that all our predictions arc in a 
good agreement with the calculations cited by Refs. 
[23-27], but our method of calculation is quite differ- 
ent from others. However, it should be also pointed 
out that both the condensates and virtualities de- 
pend on renormalizatin point /x^. We discuss this 
dependence and estimated error of in the follow- 
ing. 

A. Dependence of A^ on fj? 

As it has been mentioned in Sect. I, in QCD 
by condensates we mean the vacuum mean values 
(0 I Oi I 0) of the local operators Oj, which arise due 
to non-perturbative effects. When determining vac- 
uum condensates one implies the averaging only over 
non-perturbative fluctuations. If for some operators 
Oi the non-zero vacuum mean value appears also 
in perturbation theory, it should not be taken into 
account in the determination of the condensates. In 
other words, when determining condensates the per- 
turbation vacuum mean values should be subtracted 
in calculation of the vacuum averages. 

Separation of perturbation and non-perturbation 
contribution into vacuum mean values has some ar- 
bitrariness. Usually, this arbitrariness is avoided by 
introducing some renormalization point /i^. Inte- 
gration over momenta of virtual quarks and gluons 
in the region below /i^ is referred to condensates, 
above /x^ is referred to perturbation theory. In such a 
formulation condensates depend on the renormaliza- 
tion point yu^, (0 I Oi I 0)^2. Therefore, A^ depends 
on the renormalization point /i^, which we choose 
as 10 GeV^ as explained in the paragraph following 
Eqs(15,16). This is consistent with our conclusions 
regarding our current calculations. 

B. Estimate of Errors in A^ 

Our values of A^ for the u,d and s quarks are 
given by the quark condensate, (0 |: g(0)g(0) :| 



0), and the quark-gluon mixed vacuum condensate, 

(0 |: q[igsG^'^(7fj^[y]q :| 0). The quark condensate is 
known to about 10 %, but as discussed in Sect. II, 
approximations and models are needed to estimate 
the quark-gluon mixed condensate, and there can be 
some errors, depending on the model used. 

Although it is difficult to give an accurate estimate 
of the errors, from the results for the condensates 
and virtualities given in Tables 2, 3, and 4 one sees 
that the virtuality is estimated to about a factor of 
two. Since the parameters for Sets 1, and 2 given 
in Table 1 are more reasonable than Set 3, one may 
estimate from Table 4 that our final results for the 
virtuality are within about 20 %. 



V. CONCLUDING REMARKS 

We study the quark virtuality in the QCD vacuum 
based on the fully dressed confined quark propaga- 
tor described by DSEs. The quark virtuality is de- 
termined by the ratio of the local quark-gluon mixed 
vacuum condensates (0 |: g[«S'sG^i/iT^i/]g :| 0) to lo- 
cal quark vacuum condensates (0 |: ^(7 :| 0). The lo- 
cal quark vacuum condensate and local quark gluon 
mixed vacuum condensate are obtained by solving 
the Dyson-Schwinger Equations in the "rainbow" 
approximation with an effective gluon propagator in 
Euclidean space and the Feynman gauge. The ef- 
fective gluon propagator consists of two terms with 
two parameters: the strength of interaction x and 
its range A. The first term of the gluon propaga- 
tor simulates the infrared enhancement and confine- 
ment, and the second term matches to the leading 
log renormalization group results. Our calculated re- 
sults of local quark vacuum condensates and quark 
gluon mixed vacuum condensate are in good agree- 
ment with other theoretical model predictions such 
as QCD sum rules[24.25]^ Lattice QCD^e) and instan- 
ton modelP^l . 

Using the numerical results of our present calcula- 
tions of local quark and quark-gluon mixed vacuum 
condensates, the virtualities A^ for light quarks ( u, 
d and s ) are obtained for three different sets of pa- 
rameters X and A. The results are given in table 4. 
We find numerically that the contribution from the 
second term of gluon effective propagator in g^D'^'l, 
Eq. (12), can be neglected. The dominant contribu- 
tion to quark virtuality is from the first term of Eq. 
(12). 

In conclusion, we predict the quark virtuality us- 
ing a different method: solving DSEs, and using 
its numerical solutions Af and Bf. Our theoret- 
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ical results are consistent with all calculations of 
QCD sum rules, Lattice QCD and instanton mod- 
els. However, it should be noticed that the predic- 
tions depend on renormalization point /i^, the sep- 
aration between perturbative and non-perturbative 
part of QCD, since the vacuum condensate average 
only over non-perturbative vacuum fluctuations, and 
the perturbative contribution must be subtracted 
from any calculations. The detailed discussion on 
fi^ dependence will be published in our forthcom- 
ing paper. We believe this study is very important 
for investigation of QCD vacuum properties, and 



has many important applications both in particle 
physics and in nuclear physics. 
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